Particle-number conservation in static-path approximation for thermal superfluid 

systems 
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By applying particle-number projection to the static-path approximation (SPA), the heat capacity 
and the breakdown of pairing correlations are investigated in the thermally excited, superfluid 
systems 172 Yb, 94 Mo, and 56 Fe. For the heavy nucleus 172 Yb, the heat capacities in both the SPA 
and the number-projected SPA (NPSPA) exhibit an <S shape; the difference between the SPA and 
NPSPA heat-capacity curves is not very large and the particle-number projection thereby enhances 
the 5 shape already seen in the SPA. The temperature at which the <5>-shape of heat capacity 
curve occurs parallels the temperature of the breakdown of pairing correlations as indicated by the 
effective pairing gap. However, for the comparatively lighter nuclei 94 Mo and 56 Fe, the SPA does 
not produce an 5-shaped heat capacity on its own; only after particle-number projection the S 
shape appears in the heat-capacity curve. For 94 Mo, we compare the NPSPA result with thermal 
odd-even mass differences, which are regarded as a direct measure of the pairing gap. 
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I. INTRODUCTION 



Pairing correlations are of special importance for ma- 
ny-fermion systems such as electrons in superconduct- 
ing metals, nucleons in the nucleus, and quarks in 
color superconductivity. The Bardeen-Cooper-SchricfTcr 
(BCS) theory [l[ of superconductivity has successfully 
described the sharp superfluid-to-normal phase transi- 
tion connected to the breakdown of pairing correlations 
for an infinite Fermi system. Such a second-order phase 
transition is characterized by a discontinuity of the sec- 
ond derivative of the partition function with respect to 
some order parameter, e.g., in the present case, a jump 
of the heat capacity at a critical temperature. For finite 
Fermi systems, however, recent theoretical approaches 
0, H, 0, [f| demonstrate that thermal and quantal fluc- 
tuations are important; they wash out the discontinuity 
of the heat capacity which is obtained in the BCS ap- 
proximation, and instead produce a continuous 5-shape 
around the critical temperature. Such an S-shape has 
been found experimentally by the Oslo group [f| 0] , and 
interpreted as a signature of the pairing phase transition. 
Another fingerprint of this transition is the local decrease 
of thermal odd-even mass differences [1, [H, [l(| extracted 
from the observed level densities of a triplet of isotopes 
with neutron number N — 1, N, and N + 1, which yields 
a temperature corresponding to the one obtained from 
the S shape of the heat capacity curve. 

In finite Fermi systems such as nuclei, the BCS theory 
alone fails to provide a good approximation of thermal 
properties. This comes because particle number is not 
a good quantum number in the BCS description of the 
superfluid phase. Indeed, a BCS treatment with rigorous 
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number projection [Til such as the variation after projec- 
tion (VAP) method [111 Il2| smooths out the sharp phase 
transition and hence gives a more realistic picture, while 
the discontinuity in the heat capacity remains in, e.g . , the 
projection after variation (PAV) method [III [H, 0, [ll| . 
Exact number projection is also essential for an accurate 
description of odd-even effects seen in the heat-capacity 
curves of small superfluid systems. It is therefore de- 
sirable to investigate how the number projection affects 
thermal properties such as the S shape in heat capacity of 
nuclei llallal . In particular, it has been recently demon- 
strated [l6[, that particle-number projection within the 
BCS theory by the PAV scheme can produce an 5-shaped 
heat capacity even in the absence of a pairing-phase tran- 
sition, i.e., when assuming a constant pairing gap at all 
temperatures due to the effect of particle-number con- 
servation on quasi-particle excitations. In general, this 
method produces both an 5-shaped heat capacity due to 
the effect of particle-number projection and a disconti- 
nuity in the heat capacity related to the pairing phase 
transition. 

Inclusion of fluctuations and correlations beyond the 
BCS theory as induced by, e.g., exact particle number 
projection, can be done starting with the path-integral 
representation of the partition function. A direct ap- 
proach to evaluate thepath integral is the shell-model 
Monte-Carlo method [jj 3, However, the compu- 
tational effort is quite large, and it cannot be applied 
to the large shell- model spaces typical for heavy nuclei. 
The static-path approximation (SPA) [ID, [It], H, EH is 
therefore a useful treatment to evaluate approximately 
the partition function in finite systems with separa- 
ble interactions. In recent theoretical approaches, also 
small-amplitude fluctuations around the static path have 
been taken into account. These fluctuations give correc- 
tions to the partition function similar to the standard 
random-phase approximation around the mean field. 
Thus, the static-path plus random-phase approximation 
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(SPA+RPA) method for interactions and temperature 
regions applied in 0, [20l.|2ll. [22| gives excellent agree- 
ment with exact results 2}. In this paper, we will per- 
form the exact number projected SPA (NPSPA), since 
the SPA results for the 5-shaped heat cap acity qualita- 
tively agree with the SPA+RPA ones 1231] . F or this, we 
employ thermo field dynamics (TFD) [24],[25[, which is a 
powerful tool for describing many-body systems at finite 
temperature. 



II. METHOD 

In this work, we consider a monopole pairing Hamilto- 
nian 
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with the time reversed states k, and the pairing force 
strength G T , where t ~ tt, v stands for protons and 
neutrons respectively. Here, e k ^ T arc the single-particle 

energies and P T is the pairing operator P T = J^ fc c k T c ktT . 
l 

By means of the number-projected SPA 0, li~3l based 
on the Hubbard-Stratonovich transformation [2q ]. the 
canonical partition function is given by 
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where i?' = H—fiN, N = J^ fe cJ,Cfc is the particle-number 
operator, and (i is the chemical potential. Furthermore, 
P/v is the exact number projection defined by 
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and we define for later Xi 



A fc = v4 + A 2 witn 4 = 
G/2. Then, the thermal energy can be calculated from 
E = —dlnZpf/dp. It is now convenient to introduce 
quasiparticles by diagonalizing the Hamiltonian of Eq. 



a k 
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where the matrix satisfies unitarity Vl^W = 1. The ma- 
trices u and v in Eq. © diagonalize the pairing term of 
Eq. 11]). They are diagonal and determined by solving 
the Hartree-Fock-Bogoliubov (HFB) equations 
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where Xk is the quasiparticle energy. 

To evaluate the partition function Zm{A) in Eq. ([3]), 
we employ the treatment of number projection using 
the TFD formalism [H,II1- The TFD is known to be a 
powerful tool for a perturbative treatment within thermal 
mean- field theory [27j | . The thermal expectation value of 
an observable is thereby expressed in terms of a vacuum 
expectation value in an enlarged space. The advantage 
of the TFD is that the thermal average of an arbitrary 
operator can be easily handled in a similar manner as 
the expectation value in the zero-temperature formal- 
ism. In TFD, the Hilbert space spanned by the quasi- 
particle operators is doubled by including a fictitious 
(tilde) operator {ak, aj.}- Then, the quasiparticle vacuum 
is doubled by the corresponding vacuum |0) defined as 
afc|0) = a k \0) = 0. Next, we introduce the temperature- 
dependent operators and vacuum by the unitary trans- 
formation 
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where the generator G is given by 

G = i^^9 k \ a\a\ - a k a k 
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Here, k is the angle of the transformation. The Fock 



space is spanned by the set of operators {a[(/3), a k (/3)} 
and the vacuum (/?)). The transformations in Eqs. ([8]) 
and (OH can also be rewritten as 
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cosh 9 k 
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With this, the thermal average of an arbitrary operator 
A is expressed as the expectation value with respect to 
the temperature-dependent vacuum |0(/3)) by 



(A) = Tr (Ae-PS) / Tr (e"^) 
= (0(/3)|i|0(/3)), 



(13) 



where H is the full Hamiltonian of the system. To deter- 
mine the angle parameter 9 k , we evaluate the thermal av- 
erage of the number operator a k a k with respect to |0(/3) ) 
as 



1 Hereafter, the index r is dropped for convenience. 
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FIG. 1: (Color online) Heat capacity as a function of temper- 
ature for 172 Yb. The dashed and solid curves denote results 
from the SPA and NPSPA, respectively, the dash-dotted line 
gives the result for the independent-particle model, i.e., for 
G = 0. Experimental data (dotted line) are taken from [(|, 
Good qualitative agreement between the experiment and the 
calculations is achieved. The temperature around the local 
maximum of the Cv curve coincides well between the exper- 
iment and the NPSPA. 



Since the left-hand of Eq. (fT4"|) should give the Fermi 
distribution function, the angle is chosen as 



fk = sinlr 9 k = 



1 



e /3A fc + 1 ■ 
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Using the definition of the thermal average from Eq. 
(fT3| . the partition function Zn(A) in Eq. ([3]) can be 
rewritten as 



Z N {A) = Tr 



P N e-^ 



Tr (e~^ A A (0(P)\P N \0(P)h ( 16 ) 



with 



Tr( e -^))=in 



-(e k -fl-X k )/T 



(17) 

The expectation value (0(/3)|P/v|0(/3)) can be evaluated 
using the TFD technique [l4| in a way similar as in the 
zero-temperature formalism [28]. For instance, by ap- 
plying the general formalism of projection [l4j], we can 
obtain the expectation value of the operator e~ tS by 
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FIG. 2: (Color online) Heat capacity as a function of tem- 
perature for 94 Mo. Solid and dashed curves as in Fig. [TJ 
Experimental data (dotted line) are taken from [29l ]. There is 
good qualitative agreement of the experimental data with the 
NPSPA. In particular, the temperatures of maximum local 
enhancement of the CV curves over a Fermi gas agree well 
between the experiment and the NPSPA. 



where S is defined as the one-body operator S = <pN in 
the quasiparticle representation 
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As mentioned above, the SPA can avoid the sharp 
phase transition, which appears in the simple BCS ap- 
proximation. To explain this, we use the saddle-point 
approximation for the integral of Eq. ([2]), where we take 
into account the measure A when the maximum of the in- 
tegrand is determined. Neglecting the number projection 
for the sake of simplicity, the effective BCS approxima- 
tion (EBCS) |30( is obtained in this way. This leads to 
an effective value Aq which is determined by 
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A = G(0(/3)|P|0(/3)) + -GT/A 



(24) 
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FIG. 3: (Color online) Heat capacity as a function of temper- 
ature for 56 Fe. Solid and dashed curves as in Fig. [T] 

together with the condition for the particle number 

N = [1 - («fc - «*) tanh^Afc/2)] . (25) 

fe>0 

These equations are nonlinear and have to be solved by 
iteration in a similar way as the usual BCS equations. 
The solution of the above equations is a smooth, non- 
vanishing function of T. The second term in Eq. (|24p 
prevents the solution from falling into the normal phase. 
Now, we introduce the effective pairing gap as 

A E (T) = G(0(/3)|P|0(/?)) = Ao - GT/2A . (26) 

Here, one can see that if one neglects the second term 
in Eq. (|24p . the EBCS equation is reduced to the usual 
BCS equation. In this sense, the second term in Eq. (|2~4"|) 
washes out the discontinuity of the heat capacity and 
avoids the vanishing pairing gap at the critical tempera- 
ture in the BCS theory. Such a behavior is similar to the 
one seen in the BCS treatment with number projection 
using the VAP scheme, where the sharp phase transition 
at the critical temperature is smoothed out. 

The expectation value (0(/3)\P N \0(/3)) in Eq. (I6J) can 
now be evaluated using the M -point formula 

(0(/3)|P w |0(/3)) = ±-j\^ N (0me-^\0(f3)) 
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_ £ e ^ N {m\e- l ^ N W)), (27) 
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where (p m = 2%rn/(M+ 1) and M is the number of single 
particle states [28[. In the case M = 1, the number 
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FIG. 4: (Color online) Effective pairing gap as a function of 
temperature for 172 Yb, 94 Mo, and 56 Fe. Solid and dashed 
curves as in Fig. [1] 



rojection is reduced to the number-parity projection 0, 



III. RESULTS AND DISCUSSION 

In this work, we consider 172 Yb, 94 Mo, and 56 Fe for nu- 
merical calculations. For these calculations, we use the 
single-particle energies given by a deformed Woods- 
Saxon potential with spin-orbit interaction [3lj ]. The 
Woods-Saxon parameters are chosen such as to approx- 
imately reproduce the experimental single-particle en- 
ergies extracted from the energy levels of 133 Sn ( 132 Sn 
core plus one neutron) for 172 Yb, where Vo = 51.0 MeV, 
a = 0.67 fm, n = 0.67, A = 22.0, and r = 1.27 fm. 
For 94 Mo and 56 Fe, we adopt the same parameters as 
ones used in Ref. [10]. The doubly-degenerate single- 
particle levels with negative energies are taken outside of 
the 132 Sn core for 172 Yb, and outside of the 40 Ca core 
for 94 Mo and 56 Fe. The pairing force strengths are cho- 
sen such as to reproduce the experimental odd-even mass 
differences at zero temperature. 

Figure [T] shows the heat capacities calculated in the 
SPA and number-projected SPA (NPSPA) for 172 Yb. 
This heat capacity is obtained by dE/dT where E is 
the thermal energy given by E = —dlnZ^/d/3. One can 
see that the heat capacity exhibits an S shape around 
T = 0.5 MeV @. The number projection decreases the 
SPA heat capacity for T < 0.35 MeV and increases it 
in the region of 0.35 MeV< T < 0.7 MeV, and hence 
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enhances the S shape. It is important to note that an S- 
shaped heat capacity is already obtained within the SPA, 
without any number projection. To emphasize this point, 
we also show in Fig. [1] the heat capacity obtained from 
an independent-particle model without pairing, i.e., for 
G = 0. As it should for a non-interacting Fermi gas, this 
heat capacity depends almost linearly on temperature. 
The difference between the heat-capacity curves of the 
independent-particle model and the full model treated 
within the SPA is large and qualitative. On the other 
hand, the difference between the SPA and NPSPA results 
is relatively smaller and quantitative in nature. When 
compared to experiment, the NPSPA can describe well 
the <S-shape of the data, while the calculated C v curve de- 
viates from the experimental one for temperatures below 
0.5 MeV. This difference may be a drawback of the SPA 
approximation and of the simple model used with only a 
monopole pairing interaction and without any other in- 
teraction. As shown in Fig. [3 for 94 Mo and in Fig. [3] for 
56 Fe, however, the SPA does not produce an 5-shaped 
heat capacity by itself in those cases; the S shape ap- 
pears only in the NPSPA. A recent analysis of poles 
in the complex temperature plane [32| suggests that the 
pairing phase transition exists for mass A > 100 but not 
for A < 100. This interpretation seems to be consistent 
with the results obtained here. 

The S shape has also been discussed to be correlated 
with the breaking of nucleon Cooper pairs [3, Q . There- 
fore, we further investigate the neutron pairing proper- 
ties in the calculations 2 . In Fig. |4j we show the neu- 
tron effective pairing gap Ag in the SPA and Ag = 
G{0{(3)\P N P\0((3)) in the NPSPA, relative to their val- 
ues at T = 0. The suppression of Ag is well correlated 
in temperature with the presence of the S shape of the 
heat capacity in Fig. [IJ consistent with the results of 
Refs. [4, 9]. Thus, the <S shape in Figs. [[J d and [3] can 
be understood in terms of the suppression of the effec- 
tive pairing gap and the effects of number projection. 
Previously, we have identified the inflection point of the 
effective pairing-gap curve as the temperature at which 
the pairing transition takes place Q. As seen from the 
respective curves for 172 Yb in Fig. 0] this inflection point 
is close to the temperature 0.5 MeV at which the heat- 
capacity curves peak 3 . For 94 Mo and 56 Fe, the respective 
inflection points of 0.65 MeV and 0.9 MeV are also close 
to the local maxima of their Cv curves. 

As the experimental counterpart of the effective pair- 
ing gap, we have proposed in our previous work thermal 
odd-even mass differences (TOEMD) as a direct measure 
of the size of pairing correlations [8j] , and we have used 
them as indicators in our study of pairing phase transi- 



2 The neutron pairing properties are very similar to the ones of 
proton pairing, therefore, proton pairing is not discussed sepa- 
rately in this work. 

3 To obtain a precise estimate of the inflection point, we differenti- 
ate the effective pairing-gap curves with respect to temperature. 
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FIG. 5: (Color online) Comparison between the experimental 
TOEMD (solid line) extracted according to Eq. ([28} and the 
effective neutron pairing gap (dashed line) as a function of 
temperature for 94 Mo. 



tions in 184 W and 94 " 97 Mo [13]. We obtained in the 
one case a drastic in the other case a gradual decrease of 
the TOEMD, and we found that this signal is well corre- 
lated with the S shape of the heat capacity. According 
to Ref. [§] , the sudden decrease of the thermal odd-even 
mass differences is interpreted as a rapid breaking of nu- 
cleon Cooper pairs. Figure [5] shows the comparison be- 
tween the effective pairing gap and the TOEMD defined 



A&(Z,N,T) 



-1) 



/V 



[B t (Z,N+l,T) 



-2B t (Z,N,T) + B t (Z,N - 1,T)], 



(28) 



where the thermal energy B t is defined by B t (Z, N, T) = 
E(Z, N, T) + B(Z, N), B(N, Z) is the binding energy at 
zero temperature, and E(Z, N, T) is evaluated from ex- 
perimental level densities [3, Gil ■ The agreement between 
theory and experiment is satisfying. 

It is now interesting to discuss the significance of the S 
shape of the heat capacity. For this reason, we would like 
to recall all the available facts. In the present work, the 
S shape is obtained for A > 100 already within the SPA, 
while for lighter nuclei, to reproduce the S shape, the 
NPSPA (corresponding to a VAP scheme) is needed. In- 
terestingly, in the work of Esashika et al. [161 ] , an S shape 
could also be obtained (within a PAV scheme) when keep- 
ing the pairing gap A artificially constant. The latter two 
observations seem to indicate that the presence of the 
S shape is not necessarily related to the pairing phase 
transition, and might be connected more to the particle- 
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number projection. However, one should remember that 
in calculations without a distinct pairing force G, no S- 
shaped heat capacity has been observed. Moreover, the 
temperature for a potential pairing phase transition de- 
rived from the 5-shaped heat capacity agrees very well 
with the temperature where pairing correlations are be- 
ing suppressed (as indicated by the effective pairing gap 
Q ) , such that a simple coincidence of the two phenom- 
ena can likely be ruled out. At present, we leave the 
question of the significance of the 5-shaped heat capac- 
ity somewhat open, while we would like to point out that 
the NPSPA is certainly an important tool to investigate 
this problem, since it enables us to obtain an 5-shaped 
heat capacity for lighter nuclei in the first place. 

IV. CONCLUSION 

In conclusion, we have investigated effects of particle- 
number conservation in the SPA for 172 Yb, 94 Mo, and 
56 Fe. The particle-number projection affects the S shape 
of the heat capacity in all of these nuclei. The <S shape 
in the heat capacity of 172 Yb is produced by a coopera- 



tion of the quantum effects in the SPA with the effects 
of particle-number projection. For 94 Mo and 56 Fe, how- 
ever, the S shape of the heat capacity appears only in 
the calculation with particle-number projection, but not 
in the SPA alone. This observation is consistent with 
nuclear size effects on the existence of a pairing transi- 
tion. The 5-shaped heat capacity from NPSPA calcula- 
tions correlates well in temperature with the reduction 
of the effective pairing gap. The effective neutron pair- 
ing gap in 94 Mo is in good agreement with experimental 
thermal odd-even mass differences. Our treatment of the 
particle-number projection can be applied together with 
the angular-momentum projection in order to study the 
spin distribution of nuclear levels in such a formalism as 
suggested in Calculations are now in progress. 
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